Mirror symmetry breaking and restoration: the role of noise and chiral bias 
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The nonequilibrium effective potential is computed for the Frank model of spontaneous mirror 
symmetry breaking (SMSB) in chemistry in which external noise is introduced to account for random 
environmental effects. When these fluctuations exceed a critical magnitude, mirror symmetry is 
restored. The competition between ambient noise and the chiral bias due to physical fields and 
polarized radiation can be explored with this potential. 
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PACS numbers: 05.40.Ca, ll.30.Qc, 87.15.B- 

Enantiomers are molecules that are nonsuperimposible 
complete mirror images of each other. A remarkable fea- 
ture of Nature is that this mirror or chiral symmetry is 
broken in all biological systems, where processes crucial 
for life such as replication, imply chiral supramolecular 
structures, sharing the same chiral sign (homochirality) 
for all present living systems. These chiral structures 
are proteins, composed by aminoacids almost exclusively 
found as the left-handed enantiomers (L), and DNA, 
RNA polymers and sugars with chiral building blocks 
composed by right-handed (D) monocarbohydrates, and 
chiral amphiphiles forming membranes. This fact has 
led to the widespread perception that the presence of 
handed or chiral molecules is a unique signature of living 
systems. The emergence of this biological homochirality 
in the chemical evolution from prebiotic to living sys- 
tems is a tantalizing enigma in the origin of life, as is 
the robustness of homochirality in actual living systems, 
and is a fascinating subject that has intrigued scientists 
from diverse backgrounds. Current reviews of the origin 
of homochirality can be found in P, 0, H, 0, H, @| • Pre- 
vious hypotheses suggesting that homochirality emerged 
after the development of the primeval biological system 
, are being replaced by the widespread conviction that 
enantiomerically pure compounds are a prerequisite for 
the evolution of living species and that mirror symmetry 
breaking must have taken place before the emergence of 
life [1, HliOl- We adopt the latter viewpoint here. 

Frank introduced a paradigmatic model for sponta- 
neous mirror symmetry breaking (SMSB) and autocat- 
alytic amplification in 1953 [111]. A variant [12] of 
this open-flow reaction scheme involves the two enan- 
tiomers L and D and an achiral reactant A (kept at 
constant concentration) and the following reaction steps, 
where the k±i denote the forward/reverse (±) rate con- 
stants: Production of chiral compound (fci,fc_i): A <-> 
L, A <-!■ D, autocatalytic amplification {k2,k^2)- L + A ^ 
L + L, D + A ^ D-l-D, and mutual inhibition (k^): 
L + D ^ LD. The heterodimer LD is removed from 
the system. Frank's model contains the fundamental 
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ingredients believed to be essential for mirror symme- 
try breaking and subsequent chiral amplification [l3|. It 
can be elaborated by adding in polymerization side re- 
actions [H, m, [H, [13, [11] that can yield homochiral- 
ity in populations of oligomers. The corresponding rate 
equations expressed in terms of the enantiomeric excess 
(ee) 7] = {[L] - [D])/{[L] + [D]), the order parameter 
for mirror symmetry breaking, and the net chiral matter 
X=[L] + [D] are: 

dr]/dt = -2fciAr7/x+(fc3-fc-2)x??(l-?7')/2 (1) 
dx/dt = 2kiA - x^[k-2 + (ks - fc-2)(l - ??')/2] (2) 
+ {k2A-k^i)x. 

These rate equations are deterministic, but more realis- 
tic treatments should take noise phenomena into account. 
The nature of such fluctuations can be internal as well 
as external to the chemical system. Intrinsic statisti- 
cal fluctuations in rj about the ideal racemic composition 
[L] = [D] p^, as well as diffusion- limited noise present 
in spatially extended systems [20], are sufficient to tip 
the system over into one of its equally likely stable chi- 
ral states when fca > k-2- For prebiotic scenarios, the 
coupling of reaction schemes such as this one to envi- 
ronmental effects (e.g., meteor impacts) is crucial for de- 
termining the role of early planetary environments and 
external disruptions on the emergence, if any, of homochi- 
rality PI. 

This Letter has a two-fold purpose. On the one hand, 
we aim to establish analytically the impact of both envi- 
ronmental disturbances and chiral bias on chemical sys- 
tems that lead to SMSB. These external effects can be 
modeled stochastically and lead one to consider stochas- 
tic differential equations [11] . Recently, we developed an 
analytic perturbation method for calculating potentials 
associated with a wide class of stochastic partial differ- 
ential equations [2l|. The potential is ideally suited for 
treating symmetry breaking phenomena in nonequilib- 
rium systems. Hence, the second goal of this Letter is to 
demonstrate the computational utility of that method for 
a fundamental model of mirror symmetry breaking. The 
basic result is that ambient noise tends to restore mirror 
symmetry and homochirality is diminished, confirming 
independent numerical results [l8| . 

Tree potential. For constant A, introduce dimension- 
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less time t = {k2A — k-i)t, and we verify that when 
the rate of autocatalytic ampHfication exceeds the rate 
of monomer decay, x changes more rapidly than the 
enantiomeric excess 77. The system rapidly reaches a 
quasisteady state for x [dx/dt « 0) and then the slow 
variable 77 evolves and the full system reaches its true 
steady state [l^l- For this adiabatic regime, we then 
put X ^ X* in Eq. ([T|), where x* denotes the quasis- 
teady value for x- We define the potential V{rf) 15] via 
^ ^F{t]) = -V'{r]), and so obtain 



+ ir- 77)??^ +Vo, 



(3) 



where vq is an integration constant, and where b — 
i(fc3 - fc-2)x* > 0, r = a/b and a = hA/x* > 0. For 
the scaled potential, r is the only free variable. This 
is plotted in Fig. [T] as a function of — 1 < 77 < 1 and 
for < r < i. The absolute minima correspond to the 
asymptotic stable states of the chemical system and are 
located at 77 = ±-\/l — 2r. By varying r, we see how di- 
rect monomer production (fci > 0) tends to racemize the 
system, as the two chiral minima move continuously to- 
wards zero and coalesce at the origin when kiA increases. 
Strict homochirality = 1 holds only for kiA — 0, oth- 
erwise, kiA > implies \ri\ < 1. For r > ^, the chiral 
symmetric state ry = is the only stable solution. Gleiser 




FIG. 1: The Frank model potential V{r])/b, Eq.©, displaying 
the racemizing tendency of the direct chiral monomer produc- 
tion A — + L, A — + D. Sequence of curves from bottom to top: 
r = 0,0.2,0.4 and r = 0.5. 



an effective^ noise-corrected force -F4, as follows : 
fir/ 



dt 



(5) 



where to one-loop order in the noise amplitude A, Fj^ is 
given by [Hj 



FA{ri) = F(i^) + \^^[^^/\F{n¥+W)F^) 



VWW)+0{A^), 



(6) 



and is the real part. Thus for example, if a large me- 
teor impacts near a well-mixed prebiotic puddle or small 
pond, the bulk pond is "shaken" as a whole and a time 
dependent noise ^ {t) could provide a satisfactory descrip- 
tion of the disturbance. 

Eq.® for r = imphes F(?7) = 6ry(l - jf). The ex- 
pression under the first square root in Eq.® is [F'{rf)Y + 
F(ri)F"{ri) = b'^{l - Urj^ + 15rf). This is negative on 
the open intervals (-0.84,-0.31) and (0.31,0.84), zero 
on their endpoints, and is strictly positive elsewhere [23| . 

The one-loop effective potential is therefore given by 



Va{v) 



(7) 



where vi is an integration constant. We define 
drj 



I, = 



V[F'{v)?+F{v)F"{v). (8) 



This integral can be worked out in closed form and yields 
21, ^ _i,|V^^iV±2,^2^^,4Vi?+18(.^-l) + 8, 



7?^ - 1 



151n|2\/l5^ + 30?72 - 12| 



(9) 



valid whenever R = 1 — 127/^-f ISry^ > 0. Otherwise, from 
5R in Eq.® we have Ti = 0. Next, define T2 as follows: 



I2 = 



F(^ 



(10) 



and Walker [l8| obtained a potential qualitatively similar 
to Fig. m for a reduced polymerization model with direct 
production of monomers, which also clearly exhibits the 
racemizing tendency of such autogenic terms (see their 
Fig. la). 

The effective potential. Following Gleiser and cowork- 
ers [l^l , we couple the system to an external noise source 
^ to model random environmental effects . Applying the 
methods developed in [2i|, the corresponding stochastic 
differential equation for 77 

f = Firj) + m, mm) =ASit- f), (4) 

can be written as an ordinary differential equation with 



Since the function yJ[F'{'q)Y ^ b\{l - ?,Tf)\ then J2 

iln|7?2|+ln|l-r;2|+ci, (--L < ^ < ^) 

(11) 

-iln|7?2|-ln|l-772|+C2, (^<--i=&-i=<r;). 

Matching up at 77^ = ^ ensures continuity in I2 ■ Without 
loss of generality, we take C2 = 0. Then ci = — In i — 
21n| ~ 1.91. 

The effective potential Eq.© can be written in terms 
of these two integrals as follows: 

Va{i) = V^(r/) I2} + 0{A^), (12) 



3 



up to constants of integration used to match up the one- 
loop corrections to insure continuity. For domains over 
which R > 0, namely rj < —0.84, and —0.31 < < 
0.31 and rj > 0.84, then Ii is given by Eq.®, oth- 
erwise when i? < 0, then Ti = 0. Thus, for those 
regions over which i? < 0, the one-loop correction in 
Eg. p^ is equal to +^2^2- On the two outer intervals 
(—1, —0.84) and (0.84, 1), the one-loop correction is given 
by -4 {Ii - 12}- From Eqs.(IlII|) we can calculate this 
quantity valid on these intervals, and find that 



{Ii - I2] = -xln| — 



1 

2""' 772 ' ' 2 

+ In I 4/R+ 18(772 - 1) + 8 I 

In 12^15^ + 307/2 _ -^2| 



-ln|77 



(13) 



Whereas for the central interval (—0.31,0.31), we calcu- 
late 

1 



{X1-X2} = In |2\/;R- 12772 + 21- In |l-?72|_ci 



In 



4/R+ 18(772 - 1) 
772 - 1 



In |2V]Xr + 307/2 _ -^2| 



(14) 



Next write Va = V + {A/2)/W, then the form of the 
pure one- loop correction lS.V{rj) is completely specified 
as follows: 

SVouM + VI : (-1,-0.84) 
M^) + V2: (-0.84,-0.31) 
Ay (7,) = { ■■ (-0.31,0.31) 

l2{r]) + V2: (0.31,0.84) 
SVouM + Vi: (0.84,1). 



Here, -SVout is given by Eq.((T31), -SV^n by Eq.((Tl|) and 
I2 by Eg. pTj) . Matching up at the endpoints of the above 
intervals fixes the constants V2 = fi + 3.182, 773 = vi — 
0.001, where vi is an overall integration constant we are 
free to choose; see Eq.©. We take vi — <5Vi„(0). 

Racemization. We investigate the role that weak ex- 
ternal noise has on mirror symmetry breaking using the 
effective potential. We first scale out by the factor b, 
and evaluate Vj\,/b while varying the dimensionless noise 
amplitude < ^ <C 1. The absolute minima of the ef- 
fective potential correspond to the possible stable final 
chemical states. From the sequence of curves in Fig. [21 
corresponding to ^ = 0.0,0.05,0.1,0.2 and 0.3, we see 
that increasing the noise amplitude tends to racemize 
the system. The homochiral states I77I = 1 exist only in 
the absence of noise (bottom curve). For low levels of 
noise, the system has stable chiral states corresponding 
to 1 77 1 < 1. For noise above a critical value, the only sta- 
ble final state is the racemic solution 77 = (top curve). 
Applying a linear stability analysis to Eq.® we calcu- 
late the critical noise level Ac above which the racemic 




FIG. 2: The one-loop effective potential for the Frank 

model. The final enantiomeric excess rj, corresponds to the 
absolute minima of the potential, and decreases in absolute 
value below unity as the noise strength increases. The curves 
from bottom to top correspond to A/2b = 0.0,0.05,0.1,0.2 
and 0.3. 



state is the globally stable solution: ^ = i ^ 0.33, 
which is borne out by inspection of the curves in Fig. 
[21 New relative maxima begin to form and persist for 
A/ 2b > Ac/2b thus leading to a pair of metastable chi- 
ral states (Fig. [2l). Numerical simulations, in two and 
three dimensions IS], indicate however that the ee goes 
to zero continuously as the noise increases from zero and 
becomes strong, so these extra maxima are most likely ar- 
tifacts of the lowest order calculation. Using the nominal 
values fca « 102Ms-\ fc_2 ~ 10~^Ms-i and x* = 1 M, 
then 2b — lOOs"^, and external noises with A ^ 33s~"'^ 
would be perturbatively valid. 

Chiral bias. External magnetic, electric, gravitational 
fields, and vortex motion, as well as polarized radiation, 
can induce mirror symmetry breaking (2^ . Chiral bias 
can be studied via the potential by assigning chiral spe- 
cific reaction rates to the monomer production and auto- 
catalysis steps thus replacing ki by fcf = ki{l + ^e) and 
^ =^i(l - ^e) for i = ±1, ±2 where h = {kf + fcf )/2 
For example, e — ^ 



la [251. For example, e = ^ « 10^" for parity 
violation in the electroweak inlieractions at room tem- 
perature, where iS.E is the energy difference between 
the two enantiomers [1^. In the presence of chiral 
bias, the tree-level potential is given by (for r = 0) 



-vq, where e' 



This is plotted in Fig. [^for 1 > e' > 0: Due to the tilt, 
there are no longer racemic solutions for any bias e > 0, 
only chiral states are possible, and only one of these two 
chiral states will be an absolute minimum; see also Fig. 
lb of Gleiser and Walker . When we include noise, 
the one-loop biased effective potential is obtained by sub- 
tracting e' (7; — 77^/3) from the right hand side of Eq.(|T2]), 
and is valid up to terms of order 0{€'A) and 0{A^). At 
this lowest order, the effect of the bias is to tilt the noise 
corrected potential in the same sense as shown in Fig. 
[21 so that the sequence of noise induced minima located 
at < 77 < 1 in Fig. [21 now become the absolute min- 
ima. Due to this tilting the origin of the potential is no 
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FIG. 3: Tree potential subject to small chiral bias. Dashed 
curve corresponds to zero bias. For e' > 0, the point rj — 1 
is the absolute minimum of the potential. For the sequence 
of solid curves below dashed curve, from top to bottom, e' — 
0.1,0.2 and 0.5. 

longer locally flat (V^(0) ^ 0) for any value of the noise. 
The noise has a racemizing effect upon the biased sys- 
tem such that above a critical noise level, the effective 
potential possesses a global minimum corresponding to a 
weakly chiral state. Thus for example, we calculate that 
for e' = 0.001,0.01 and 0.1 then ^ = 0.311,0.318 and 
0.388, and the corresponding enantiomeric excesses are 
rj w 0.10,0.15 and 0.20, respectively. For increasing chi- 
ral bias, ever stronger noise levels are "tolerated" before 
homochirality |r/| = 1 is erased. A detailed account of 
noise and chiral bias on mirror symmetry breaking will 



be provided elsewhere. 

In this Letter we applied the stochastic field theory 
formalism of [2l| to study the emergence of chirality in 
a key model of SMSB in chemistry in which environ- 
mental effects are modeled by external noise. We fo- 
cused on the Frank model due to the central role it plays 
in theoretical app roaches to mirror symmetry breaking 
[i [III, [11, H, [S [M [B H [SO, III, [il 111 ■ By strictly 
analytic means we verified that weak noise racemizes the 
system, erasing homochirality. This is a perturbatively 
valid key result, confirming the previous numerical re- 
sults obtained by Gleiser and coworkers [l^. We also 
studied the competition between chiral bias and external 
noise and verified that stronger noise levels are required 
to racemize the system in the presence of bias. We as- 
sumed well-mixed conditions (zero dimensional systems) , 
but the analytic method ^] enjoys the flexibility to in- 
clude diffusion in d-dimensions and spatially dependent 
noise terms. A preliminary study of the d = 2 potential 
indicates that the results presented here carry over when 
spatial dependence is included. The important role that 
fluctuation phenomena, noise and chiral bias play in the 
origin of homochirality can therefore be analyzed in an 
elegant and systematic way. 
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